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ABSTRACT 



Test particles infalling from infinity onto a compact spherical star with a mildly super-Eddington luminosity at its 
surface are typically trapped on the "Eddington Capture Sphere" and do not reach the surface of the star. The 
presence of a sphere on which radiation pressure balances gravity for static particles was first discovered some twenty 
five years ago. Subsequently, it was shown to be a capture sphere for particles in radial motion, and more recently 
also for particles in non-radial motion, in which the Poynting-Robertson radiation drag efficiently removes the orbital 
angular momentum of the particles, reducing it to zero. Here we develop this idea further, showing that "levitation" 
on the Eddington sphere (above the stellar surface) is a state of stable equilibrium, and discuss its implications for 
Hoyle-Lyttleton accretion onto a luminous star. When the Eddington sphere is present, the cross-section of a compact 
star for actual accretion is typically less than the geometrical cross-section ttR 2 , direct infall onto the stellar surface 
only being possible for relativistic particles, with the required minimum particle velocity at infinity typically ~ 1/2 the 
speed of light. We further show that particles on typical trajectories in the vicinity of the stellar surface will also be 
trapped on the Eddington Capture Sphere. 
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1. Introduction 

In Newtonian theory both gravity and radiative flux in 
vacuo diminish like 1/r 2 with the distance from the cen- 
■ ter of the star. The luminosity in units of the Eddington 
S-h \ luminosity, L{r) / L^dd, at any distance from the star is 
^ equal to its value at the surface, L(R)/LBdd, and is 
therefore everywhere either sub-Eddington, Eddington, or 
super-Eddington. However, in Einstein's general relativ- 
ity the radiative force diminishes more strongly with dis- 
tance than the gravitational force, the rcdshiftcd luminos- 
ity L(r) decreasing as 1/(1 — 2Ro/r), and a static bal- 
ance of forces is achieved at that radius r at which L(r) = 
-^Edd(l — ^Rg/ 7 ")^ 1 ^ 2 ■ For this reason, radiation may be 
super-Eddington close to the star, but sub-Eddington fur- 
ther away, reaching the Eddington value at the "Eddington 
Capture Sphere" , whose radius is ( |Phinney 1987[ ) 

2R G 



Rq defined in Eq. ([3]). Abramowicz, Ellis & Lanza (1990) 
(hereafter AEL) have rigorously shown by analytic calcuia- 
tions that in the combined gravitational and isotropic radi- 
ation fields of a spherical, non-rotating, compact star, radi- 
ally moving test particles (with proton mass, and Thomson 
cross-section for photon momentum absorption) are cap- 
tured on this sphere and "levitate," i.e., remain at rest. 



^Edd 



1 



1 2RgY f L{R)\ 

R J V ^Edd / 



(1) 



At this radius the gravitational and radiative forces balance 
(for a static, optically thin hydrogen shell). The formula is 
in Schwarzschild coordinates, with the gravitational radius 



Oh, Kim & Lee (2011) have shown that when the lumi- 
nosity at the surface of the star is mildly super-Eddington, 
particles with non-zero orbital angular momentum can also 
be captured at i?Edd (see also lBini et al. 2009]) . This is be- 
cause the Poynting-Robertson radiation drag acts as an ef- 
fective torque, and for a wide class of trajectories reduces 
the orbital angular momentum to zero. 

In this paper we further discuss the phenomenon of the 
Eddington Capture Sphere by considering two particular 
cases, one corresponding to initially quasi-circular motion 
close t o the star, and the second b eing rather similar to the 
classic Hoyle & Lyttleton (1939) accretion model. We also 
consider the question of stability ind speculate about possi- 
ble astrophysical manifestations of the Eddington Capture 
Sphere. 

All numerical results presented here were obtained with 
the Dormand-Princc method, which is a fourth-order ac- 
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curacy, adaptive step-size Runge-Kutta type integration 
method. 



2. Equations of motion 

We carry out all calculations in the Schwarzschild metric 
using standard Schwarzschild spherical coordinates, 



t = time, r — radius, 

8 = polar angle, 4> = azimuthal angle, 



and extend the treatment of AEL to non-radial motion. 
Spherical symmetry of the problem assures that the tra- 
jectory of a test particle is confined to a single plane for 
any set of initial conditions. Without loss of generality, we 
assume that this is the equatorial plane 9 = it/2 in the 
Schwarzschild coordinates. Note that the stress-energy ten- 
sor of radiation calculated in the stationary observer 
tetrad by AEL does not depend on the particle motion. 
Therefore, we may use here the AEL formulae for 
with no change. We use dimcnsionlcss radius x, and dimen- 
sionlcss proper time r defined by the scaling of the line 
clement ds, as well as the dimcnsionlcss radius of the star, 



dr = 



ds 



r 

Rg 



X 



R 

— — , with R G 
Kg 



GM 



(3) 



and introduce B = 1 — 2/x. 

Assuming that the test particle is capturing momentum 
at a rate proportional to the radiative flux in its rest frame, 
thus suffering a rest-frame force of (<j/c)x flux, the particle 
trajectory is described by two coupled, second order differ- 
ential equations: 



d 2 x 



dr 2 ttI(R)X 2 
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with the time component of the particle four velocity 
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and 



e = BT^iv*) 2 + B-'fWW f^Y 
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dx 



The parameter k is the surface luminosity in Eddington 
units, 



-^Edd 

where ^Edd = &rGMm p c/a . 
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Fig. 1. Equilibrium on the Eddington Capture Sphere. 
Trajectories of particles returning to the Eddington sphere 
{X = 5, k = 1.5) are shown for four initial velocity per- 
turbations: = xdip/dr = 0.12, v r = dx/dr = (thick 
green continuous line), = 0.20, v r = (thin blue contin- 
uous line), = 0.12, v r = 0.10 (thick red dashed line) and 
= 0.12, v r = -0.10 (extra thick violet dashed line). Also 
shown is a tangential escape trajectory (in thick brown) 
with initial = 0.32, v r = 0. 

3. The capture sphere 

The system of Eqs. (U]), © allows a static solution 
x(t) = x E dd for 



k= 1- 
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X 



(9) 



where x^dd = REdd/Rc, cf. Eq. ([T]). In the case of purely 
radial motion it was shown (AEL) that a critical point of 
nodal type is present at x = £Edd, resulting in a stable 
equlibrium state of captured particles. A corresponding set 
of equilibrium points can be recognized in the system of 
Eqs. ((4]), ([5]). The set consists of points with x = x^dd, 
and any value of and 9. (From Eq. [S] it follows that if 
d(j)/dt = 0, then d 2 (j)/dt 2 = 0, and Eq. [3] reduces to the ra- 
dial equation of AEL.) We refer to this set as the Eddington 
Capture Sphere. The radius of the sphere changes from 
^Edd = X to x'Edd = oo as fc increases from (1 — 2/X) 

to (i - 2/xy 1 . 



(7) 4. Stability of equilibrium 



Radial velocity perturbations were already discussed in 
AEL. In the case of general velocity perturbations, we ob- 
serve stability in numerical experiments: if the azimuthal 
component of the velocity of a particle located on the cap- 
ture sphere is slightly (or not so slightly) perturbed from 
zero value, the particle changes its location but returns to 
the stationary state somewhere else on the sphere. In this 
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sense, motionless "levitation" at a distance i?Edd from the 
star is a stable state. Such behavior is presented in Fig. [1] 
for k = 1.50, X = 5.00 (thick continuous black line), cor- 
responding to XEdd = 10.5 (thin dashed black line). We 
only show a small subset of the trajectories that we have 
computed, and for clarity we only show orbits with large 
initial velocity perturbations ~ 0.1c. For smaller perturba- 
tions the result is the same. 

In fact, one can compute the minimum escape veloci- 
ties from the Eddington Capture Sphere, these correspond 
to radial trajectories and are v r + = 0.23 for outward motion 
(escape to infinity) and v r _ = 0.51 for inward motion (cap- 
ture by the star), for the parameters considered. The mini- 
mum escape speed for trajectories tangent to the Eddington 
sphere (the second cosmic speed) is = 0.32, the corre- 
sponding trajectory is also exhibited (in brown) in Fig. [T] 
One of the reasons for these unexpectedly large values of 
the escape velocities is the drag experienced by the particle 
even in purely radial motion. 
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5. Capture from circular orbits in the stellar vicinity 

In the absence of radiation, a particle with appropriate an- 
gular momentum will follow a circular orbit around the 
star, but in a strong radiation field the particle will spi- 
ral inwards, since the radiation drag continuously removes 
angular momentum. 

From Eq. Q, the azimuthal speed, v c i TC , correspond- 
ing to a trajectory initially tangent to a circle, may be 
derived. Figure [5] shows some (non-escaping) trajectories 
ending on the Eddington sphere at XEdd = 10.5 (X = 5.00, 
k = 1.50), for particles with an initial zero value of ra- 
dial component velocity, dx/driQ) = 0. Corresponding ra- 
dial and angular velocity plots are presented in Fig. [3] The 
trajectory shown with a thin continuous red line starts at 
x '= 15.0 with = v circ = 0.00344. The thick dashed blue 
one and the thick continuous green one start at xq = 8.00 
with = 0.00344 and = 0.0300, respectively. They arc 
all (eventually) attracted to the Eddington capture sphere, 
and in the end the approach is almost radial due to the ef- 
fective removal of all angular momentum by the radiation 
field. 



6. Hoyle-Lyttleton accretion 

We consider a situation resembling the model of a star 
moving through a unifo rm gas cloud, introduced by 
Hoyle & Lyttleton (1939) In this Section we take the star 
to have radius X = 6. 

In Fig. [4j numerically calculated trajectories of in- 
coming particles are shown. The parameters used were 
k = 1.49 (which implies XEdd = 151), and the initial 
velocity at infinity Woo = 5.00 • 10~ 3 c (parallel to the x 
axis) for all particles. The initial value of rcos</> was 
taken to be — 5000i?c and the impact parameters (initial 
rsin</>) of presented trajectories were taken from the set 
{0, 30i?c, 60i?c, 1500-Rg}- Radial and angular compo- 
nents of velocity are presented in Fig. [5] for three selected 
particles, as a function of r. The reduction of velocity along 
the trajectory is attributed to the influence of radiation, 
bringing the particles to a halt at x = x^dd- Close to the 
equilibrium sphere, motion is almost radial and trajectories 
become perpendicular to the surface of the sphere. 



Fig. 2. Different trajectories for purely tangential initial 
velocity. 




t, 0.015- 




Fig. 3. Evolution of radial and angular velocity for trajec- 
tories in Fig. [2] 



For a certain range of kinetic energies, the incoming par- 
ticles can penetrate the capture sphere without reaching 
the stellar surface. Inside the sphere, radiation dominates 
over gravity and the particles will be pushed out until they 
reach the surface of the capture sphere again, and come to 
rest there. This is illustrated in Fig. [6j with the parameters 
chosen to be Voo = 0.200c, X = 6.00, k = 1.45, yield- 
ing XEdd = 30.5. The initial value of rcos0 was equal to 
— 2000i?c an d the impact parameters of presented trajec- 
tories were taken from the set {0, 2Rq, 4Rg, 50i?c}. 
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Fig. 4. The Hoyle-Lyttleton accretion problem for X=6.00, 
and k=1.49. The trajectories of particles for Voo = 5-10 _3 c, 
are shown. None of the test particles is accreted by the star, 
they are either trapped on the Eddington sphere or escape 
to infinity. 



;ure Sphere around luminous stars 
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Fig. 6. Trajectories of particles for v x = 0.2c with (top, k 
= 1.45) and without (bottom, k = 0) radiation influence. 
The stellar surface is indicated by the inner circle (thick, 
continuous line), and the Eddington sphere by the dashed 
circle. 



(impact parameters equal to 25.0 and 22.4) are drawn with 
thick dashed lines in the upper and lower panel, respec- 
tively, of Fig. [6j 

Park (2012)| in general, ra- 



As noted by Oh, Chun 
diation drag is responsible for capturing particles from a 
larger area than in the case of accretion with no radiation. 
However, for the trajectories considered so far, no parti- 
cle falls onto the surface of the star. For particles actually 
reaching the stellar surface, the converse statement is true, 
they come from a smaller cross-sectional area than in the 
case of accretion with no radiation, and in the cases con- 
sidered here, smaller even than the cross-sectional area of 
the star irR 2 . 



Fig. 5. Evolution of radial and angular velocity for three 
trajectories from Fig. 0] with impact parameter equal to 
(thick continuous green), 500 (thin continuous red) and 
1000 (dashed blue line). 



The largest initial impact parameter of the captured 
particles was approximately equal to 1400 and 25.0 in the 
two presented cases (Fig.|4]and upper panel of Fig. [6j, while 
the critical impact parameters of Hoyle-Lyttleton accretion 
in the absence of radiation would have been equal to 849 
and 22.4, respectively. Critical trajectories separating the 
trajectories of escaping particles from the captured ones 



6.1. Direct capture onto the star 

For a sufficiently low impact parameter, and a sufficiently 
high incoming velocity at infinity, a test particle may pen- 
etrate the Eddington sphere deeply enough to reach the 
stellar surface. The minimum value of Woo required to reach 
the star in the case corresponding to Fig. [4j i.e., k — 1.49 
is shown as a function of the impact parameter in the top 
panel of Fig. [7] (the monotonically growing function de- 
picted by the thick black line), for these velocities the par- 
ticles settle on the stellar surface with zero velocity. This is 
a case of accretion not leading to any release of energy on 
the stellar surface. If the velocity at infinity is higher, the 
test particle will deposit some kinetic energy at impact on 
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Fig. 7. Regions in the parameter space corresponding to 
escaping particles (region denoted by "oo" above the blue 
curve), particles ending up on the Eddington Capture 
Sphere (the shaded region between the thick black curve 
and the blue curve, denoted by "ECS"), and particles im- 
pacting the stellar surface (region below the thick black 
curve, denoted by "ss" ) , for two values of luminosity, k = 
1.49 (top), and k = 1.25 (bottom). The radius of the star 
is X=6 for all plots in this Section. The impact parameter, 
b, is in units of Rq, and in units of c. Particles with 
initial parameters on the thick black curve settle on the 
surface of the star with zero velocity. The thin black curves 
correspond to particles impacting with (from left to right): 
v r (R) = 0.05,0.10,0.15,0.20. 

the stellar surface. However, the deposited kinetic energy 
will be a tiny fraction of the kinetic energy at infinity, as 
the impact velocities are much lower than Voo (thin black 
lines in Fig. [7]). This outcome is just the opposite of what 
occurs when radiation drag and pressure may be neglected. 

Other cases (fc = 1.45, and k = 1.25) are illustrated in 
the bottom panel of Fig. [7] and in Fig. [3] The bottom panel 
of the latter figure illustrates the importance of radial radia- 
tion drag. The top curve illustrates the changing velocity of 
a test particle accelerated by gravity and decelerated by the 
radiation pressure [i.e., the term proportional to the (t)(r) 
component of the stress tensor] , but not suffering from radi- 
ation drag. The bottom curve shows the velocity computed 
with our full equations — clearly the resistance offered by 
the radiation field slows down the particle very effectively. 
Incidentally, the top curve illustrates another rclativistic ef- 
fect. For a rapidly moving particle the radiation pressure is 
enhanced by the ratio of E/(mc 2 ) and it overwhelms the 
gravitational pull, so that the test particle is decelerated 
well outside the Eddington sphere. 

7. Conclusions and discussion 

Wc have shown that a luminous star in Schwarzschild met- 
ric can capture test particles from a wide class of orbits 
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Fig. 8. The case of k=1.45. Top: Same as Fig. 7. Bottom: 
The velocity as a function of radius for radial infall (b = 0) 
and Voo = 0.82c, with and without radiation drag (bottom 
and top curves, respectively). 

onto a spherical surface of radius larger than the star itself. 
The particles come to rest on this sphere, because radia- 
tion drag eventually removes all of their angular momen- 
tum, as shown by Oh, Kim, & Lee (2011). The particles 
remain suspended above the surface of the star, because 
every point on the capture sphere corresponds to a po- 
sition of equilibrium, the equilibrium being stable in the 
radial direction (Abramowicz, Ellis, Lanza 1990) and neu- 
tral in directions tangent to the sphere. The radius of this 
"Eddington Capture Sphere" depends on the ratio of the 
stellar luminosity to the Eddington luminosity (Eq. [1]). 
Thus, accreting compact stars, as well as stars which eject 
matter, may be surrounded by a shell of matter at rest, as 
long as the luminosity of the star is close to the Eddington 
luminosity (super-Eddington on the stellar surface). 

The trajectories of particles falling from infinity in the 
gravitational field of a mildly super-Eddington compact 
star were investigated in Section 6. Wc showed that ra- 
diation drag strongly reduces both the azimuthal and the 
radial component of test particle velocities. This leads to a 
moderate enhancement of the capture cross-section by the 
Eddington sphere, relative to the cross-section for accre- 
tion by a non-radiating star, and a drastic suppression of 
the cross-section for actual accretion onto the stellar sur- 
face. Radiation drag also leads to a dramatic decrease in 
the kinetic energy deposited at the surface. 

It is interesting to speculate that the presence of the 
Eddington Capture Sphere, an effect of Einstein's general 
relativity, may play a role in the formation and ejection 
of shells of matter by massive stars in the final stages of 
their evolution, such as Wolf-Rayet stars, or the luminous 
blue variables (LBV stars). In principle, the capture sphere 
can exist at any radius outside the star, if the stellar lumi- 
nosity has the appropriate value. However, the larger the 
radius of the sphere (in units of the Schwarzschild radius), 
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the more finely tuned the luminosity has to be — for a IORq 
(5 Schwarzschild radii) capture sphere the luminosity at in- 
finity has to be 0.95-LEdd, while for a 500 000 Rg sphere 
the luminosity has to be equal to the Eddington value to 
an accuracy better than one part in a million. While our 
discussion was couched in terms of the Thomson scattering 
cross-section, similar results will hold if other scattering 
or absorption processes contribute to momentum transfer 
from the radiation field to the particles. The capture sphere 
is located at that radius at which radiation pressure on 
matter at rest balances the gravitational pull of the star. 
Outside the sphere, the radiation pressure is too weak to 
balance gravity, inside the sphere it overcomes the gravita- 
tional pull. 



It is expected that a more likely application of the 
Eddington capture sphere will be found in the interpre- 
tation of the behaviour of accreting neutron stars in low 
mass X-ray binaries. In particular, in the Z sources, the 
inner radius inferred from both the kHz QPOs and the 
accretion-disk spectral component have been determined 
to vary rapidly at a lumi nosity which is thought to be close 
to the Eddington value (fYu, van der Klis, fc Yonker, 2001| 
|Lin, R emillar d7"&: Homan, 2009[ ). In particular, the kHz 
QPOs vary with the normal branch oscillation (NBO) 
phase, which may indicate radiation induced rapid drifts of 
the innermost accretion flow. Perhaps these findings could 
be understood in terms of the expected behaviour of ac- 
creting matter accumulating on (a sector of) the Eddington 
sphere. In the non-spherical geometry of accretion disks, ra- 
diation scattered in non-radial directions can escape from 
the system. Therefore the accreting fluid can "levitate" 
above the star only as long as it is optically thin. As more 
and more matter accumulates at the capture radius of Eq. 
(1), the shell eventually becomes optically thick and the 
outer layers of the fluid outweigh the radiation pressure 
support, leading to rapid accretion of matter and to the 
consequent evacuation of the region near the Eddington 
capture sphere, allowing the process of accumulation to be- 
gin again. One would also expect the accreting fluid to be 
spread out to higher altitudes (lower and higher values of 
the polar angle 0) than those subtended by the inner parts 
of the optically thick disk, as the incoming fluid forms a 
"puddle" on the surface of the Eddington capture sphere. 
A more detailed understanding of actual behaviour of Z 
sources would require relaxing the assumption of spherical 
symmetry, on which all the calculations in this paper were 
based. Comparison of the theoretical expectations with the 
observed NBO and its related flaring branch oscillation 
(FBO) would be very interesting. 
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In the atoll sources, which are less luminous, the 
Eddington luminosity is attained during X-ray bursts. It 
would be worthwile to investigate whether the Eddington 
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